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We investigate quantum walks in multiple dimensions with different quantum coins. We augment 
the model by assuming that at each step the amplitudes of the coin state are multiplied by random 
phases. This model enables us to study in detail the role of decoherence in quantum walks and to 
investigate the quantum-to-classical transition. We also provide classical analogues of the quantum 
random walks studied. Interestingly enough, it turns out that the classical counterparts of some 
quantum random walks are classical random walks with a memory and biased coin. In addition 
random phase shifts "simplify" the dynamics (the cross interference terms of different paths vanish 
on average) and enable us to give a compact formula for the dispersion of such walks. 

PACS numbers: 03.67.-a, 05.40.Fb 



I. INTRODUCTION 

The concept of quantum walks (QW) has been intro- 
duced (see Ref. pj) in order to explore how the intrinsi- 
cally statistical character of quantum mechanics affects 
statistical properties of quantum analogues of classical 
random walks. In particular, an example of a random 
process is a Markov chain such that the position value 
x € X is iteratively updated, given by the transition 
probability P{x\y). 

Quantum walks have been studied in connection with 
novel quantum algorithms: The instances thereof were 
provided in Ref. [2( (the quantum walk algorithm on the 
hypercube with complexity 0(^/n)) and in Ref. Q (the 
quantum walk algorithm for subset finding) . The former 
uses the quantum walk on the hypercube, while the latter 
uses the quantum walk on bipartite graphs. Quantum 
walks on bipartite graphs were analyzed in Ref. 

Various aspects of QWs have been studied in detail re- 
cently (for a review on QWs see Ref. 0). In particular, 
Aharonov et al. have presented an analytic description 
of discrete quantum walks on Cayley graphs A spe- 
cial case of a Cayley graph, the line, was asymptotically 
analyzed in Ref. 7]. It has been shown that, unlike clas- 
sical random walks, the probability distribution induced 
by quantum walks is not Gaussian (with a peak around 
the origin of the walk), but has two peaks at positions 
±-^!=, where n is the number of steps. As a result the 
dispersion of probability distribution for quantum walks 
grows quadratically, compared to linear growth for clas- 
sical random walks. The role of decoherence in quantum 
walks has been analyzed by Kendon et al. Ref. [a, Oj 

Quantum walks are intrinsically deterministic pro- 
cesses (in the same sense as the Schrodinger equation is 
a deterministic equation). Their "classical randomness" 
only emerges when the process in monitored (measured) 
in one way or another. Via the measurement, one can 
regain a classical behavior for the process. For instance, 
by measuring the quantum coin, the quadratic disper- 



sion of the probability distribution reverts to a classi- 
cal, linear dispersion. If the quantum coin is measured 
at every step, then the record of the measurement out- 
comes singles out a particular classical path. By averag- 
ing over all possible measurement records, one recovers 
the usual classical behavior 0, • Instead of measur- 
ing the quantum coin after each step, an alternative way 
to regain classical randomness from a quantum walk is to 
replace this coin with a new quantum coin for each flip. 

After n steps of the walk one accumulates n coins that 
are entangled with the position of the walking particle. 
By measuring a set of n quantum coins, one could re- 
construct a unique classical trajectory and by averaging 
over all possible measurement outcomes, one once again 
recovers the classical result. 

These two approaches to regaining classical behavior 
from the quantum walk have been contrasted in a recent 
work by Brun et al. [Tl| . This comparison has been 
studied for the particular example of a discrete walk on 
the line. 

In the present paper we analyze the quantum-to- 
classical transition using random phase shifts on the coin 
register. In Sec. II we give an introduction to the quan- 
tum walk model. In Sec. Ill (part A) we augment the 
model by random phase shift dynamics and present the 
solution in terms of path integrals. It turns out that on 
average the interference of amplitudes of different paths 
is zero and we derive the formula for the dispersion of 
the mean probability distribution in compact form. Wc 
contrast the dynamics of quantum walks with two coins 
(permutation symmetric and Fourier transform) with the 
dynamics of classical random walks and find an equiva- 
lence between the two (considering the possibility that 
the CRW has memory and a biased coin). In part B 
of Sec. Ill we provide the numerical results of the prob- 
lem. In particular, we briefly analyze a situation in which 
phases of random shifts are distributed according to a 
normal distribution that is peaked around the phase zero 
and with the dispersion a . When the dispersion is zero, 
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(i.e. a = 0) we recover the QW, while for large a, we ob- 
tain a uniform distribution on the interval [— 7r, 7r] and the 
CRW is recovered. In between we can observe a contin- 
uous quantum-to-classical transition of quantum walks. 
In Sec. IV we present our conclusions. Some technical 
details of the calculations can be found in Appendix A. 



II. QW IN MULTIPLE DIMENSIONS 

Let us first define a quantum walk in d dimensions 
- i.e. on the lattice Z d . The quantum walk is gener- 
ated by a unitary operator repeatedly applied on a vector 
from a Hilbert space 7i = Hx <8> 'Hd- The Hilbert space 
Hx = span{|x) : x S Z d } is called the position Hilbert 
space. For Vx, y S Z d we define the usual scalar prod- 
uct x • y = ^2j = i X-jYj an d the norm |x| = y/x ■ x. In 
the following, the distance between the vertices is a di- 
mensionless quantity, with the distance between adjacent 
vertices equal to 1. 

There are 2d vectors e a e Z d such that |e a | = 1. 
The space Hd = span{|a) : a = 1, ...,2d} is spanned 
by states isomorphic to e„. 7Yzj is called the direction 
Hilbert space. In the following we set D = {1, . . . , 2d}. 

A single step of quantum walk is generated by the uni- 
tary operator U such that U = S(l ® C), where 



*=££i 



e a)(x| 



EE J2 T a ® P * S C 2 - 1 ) 



P a ee |a)(a| and C is any unitary operator. The oper- 
ator S changes the state of the position register in the 
direction a, while the coin operator C operates on the 
direction register. For simplicity we consider the permu- 
tation symmetric coin 



C 



r t 
t r 



(2.2) 



The quantum walk is generated by a sequence U n \'ipo}, 
where \ipo) is some initial state. For simplicity, we assume 



where 



>\D\ 



3 \s) , (2.3) 
We also assume the so- 



called Grover coin [12j , which is a specific instance of the 
permutation symmetric coin in Eq. 12. 2f) . described by 
the operator 



C G = 2| S )(s|-l 



(2.4) 



In order to find the eigensystem of U , we switch to the 
translationally symmetric basis 0- We set 



£ e* x |x) 



(2.5) 



where k S K d . By virtue of the inverse Fourier transform 
we obtain 



1 



(2ir) d 



c |^k)d d k , 



(2.6) 



with 1 0k) are eigenvectors of the translation operator in 
the a-th direction, i.e. 



e- lk - e "|0k) 



(2.7) 



where T a = J2xez d |x+e a )(x|. By applying the evolution 
operator, we obtain 

t/|<M® |x) = s\4> k ) ® c\ x ) 

= 5>-* k - e »|0 k )®|a>(a|C|x> (2g) 

a£D 

= \4>k) ® A fe C|x) , 



where 



A fe = E e " k ' e i a >( a i 



(2.9) 
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In order to simplify the notation in what follows we will 
denote projectors \a)(a\ as P a - To find the eigensystem of 
U, we need to find the eigensystem of AfcC Equivalcntly, 
we need to evaluate (AkC) n \x)- 

We first use the Grover matrix Cq in Eq. (|2.4|l . In 
order to find the power of the matrix (AkCc)™ we prove 
the following lemma: 

Lemma 1. Let T> = {\a}} be the orthonormal basis 
of a Hilbert space and Cq — 2P S ~ 1, where \s) — 

-/pEaGDla)' Then 



with p n 



<\D\ 



{P a C G ) n = Pn \a)(s\ + q n P a , 

(m ~ x ) and q ' 1 = ~ 



2 - 1 

\D\ 1 



Proof. We denote P a C G = -/=|a)(s| - Pa = Po\a)(s\ + 

q P a - Setting (P a C G ) k = Pk\a)(s\ + q k P a we get that 
Pk+i = -Pfc(-r^T — l) and(7 fe+ i = q k (jj^-l). By induction 
we immediately obtain the result. □ 

From Eq. (|2.9(l we see that with the Grover coin 

(A k c G r = (j2 



aeD 



= e-*-^+-+ e ^P ai C G ...P an C G . 

(oi,...,njeD" 



(2.10) 



By induction, the expression P ai C G . . . P an C G from 
Eq. I|2.10|l can be rewritten as 
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Lemma 2. 

P ai C G . . . P an C G 



(-1)* 



|£)|(2n-l)/2 



n-l 



x II(l' D l*«i.-i+ I - 2KV\D\\ai)K\ - 2| ai )( S |) . 

(2.11) 

The product in Eq. i|2.11|l is taken to be 1 for n = 1. 

Alternatively, in Eq. I|2.1()(l , the last line can be rewrit- 
ten as 



P ai Cc ■ ■ ■ Pa n Cc 

= n (Pa^c G ) m >...(p a(k) c G r* 



(2.12) 



mH hwifc — ?! 



where € D. According to Lemma ^ all the terms 
in the product in Eq. (|2.12|l can be expressed as the 
linear combination of \a^){s\,P a (j). Since for j ^= f => 
{a^'\a^ >) = 0, we get the result 



(A k C G )"= J2 e- ik < e aW+-+ e aW 

partition 

x ( S(mi, ■ ■ .,m k )\ai)(s\ + T(m 1 , . . . 



> m k )Pa n 



(2.13) 



The expressions for S and T are given by relations 



S (mi, 
T(mi, 



,m k ) 
,m k ) 



Pmi 
Pmi 



■Pmk-iQmk , 



(2.14) 
(2.15) 



where we use the notation of Lemma ^ The coefficients 
mi, . . . ,m k give the partitioning of the integer n such 
that mi + • ■ ■ + m k — n. The "partition" in Eq. (|2.13|) 
means the summation over all such partitions. 

Starting with the initial state \ipo) — |0)<£>|s) and using 
the expression H2.11fl . we obtain 



(AkC G ) n |s) 



— ik-(e ai H he 



£ 

■-,0*71, 

L 

[n^iwx- 2 



71+ 1 



(ai,...,a„)6_D 
n-l 



|£)|(2n-l)/2 



(2.16) 



3=1 



This equation takes the singular form for \D\ =2 (i.e. , 
for one-dimensional quantum walk on line) such that only 
the summands for which all elements {ai,...,a n } are 
distinct, contribute to the total sum. As a consequence, 
the sum in Eq. (|2.1(il) when \D\ = 2 is zero for n > 2. 
But this case is special in that the coefficient r of the 
Grover matrix is zero. From now on we will consider the 
dimension of the lattice to be equal or larger than 2, so 
that, \D\ > 4. 



The expression (12.16(1 is symmetric with respect to the 
permutation of elements, i.e. , (aj\(A^Co) n \s) has the 
same value for any a,j € D, n G Z. A value of the right- 
hand side of Eq. 12.16fl depends on the term 

n-l 

S (oi,...,On)= l[(\D\6 aj ,a j + 1 -2). (2.17) 

i=i 

Obviously, |So(<Xi, . . . , a n )\ is maximal for a\ = ■■■ = 
a n . More precisely, if aj = Oj+i for j = 0, . . . , n — 1 in 
Eq. l|2~T7jl . then 



|3 (a 1 ,...,a„)|=O((|£>|-2) fc 2 



k c\n—k\ 



(2.18) 



and there are O (\D\ n ) such terms in the sum of 
Eq. (gHHJ}. Now Eq. (f^TTBf) takes the form 

(A k C G )» = 

= £ e-^^-^E Q (ai,...,a n )\ai) 

1 1 (o 1 ,...,o B )£D» 

(2.19) 

In what follows we will compare the quantum walk de- 
scribed by Eq. l|2.8|l with the quantum walk with random 
phase shifts. 



III. QW WITH RANDOM PHASE SHIFTS 
A. Analytic Results 

Quantum walks differ from classical random walks 
(CRW) in many respects. One of the main differences 
is that dispersions of probability distributions of CRW 
grow linearly with the number of steps while for QW the 
dispersions grow quadratically 0- In what follows we 
will show that introducing random phase shifts (RPS) at 
each step of the evolution causes the QW to behave more 
like a classical random walk. The reduction of the QW 
to the CRW has been discussed in Refs. 0, [n| . The au- 
thors of these papers have discussed two possible routes 
to classical behavior for the discrete QW on a line. First, 
the QW-to-CRW transition has been considered as a re- 
sult of decoherence in the quantum "coin" which drives 
the walk. Second, higher-dimensional coins have been 
used to "dilute" the effects of quantum interference. The 
position variance has been used as an indicator of classi- 
cal behavior. It has been shown that the multi-coin walk 
retains the "quantum" quadratic growth of the variance 
except in the limit of a new coin for every step, while the 
walk with decoherence exhibits "classical" linear growth 
of the variance even for weak decoherence. 

In what follows we will utilize a different approach to 
analyze the QW-to-CRW transition. In Ref. [ll| the au- 
thors used a CP-map on the coin degree of freedom to 
simulate the effects of decoherence on the quantum walks 
in 1 dimension. If the CP map is pure dephasing, then 
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the dispersion of the probability distribution is asymp- 
totically linear. Our approach is new in two respects: 
Firstly, we generalize the problem to an arbitrary num- 
ber of dimensions; secondly we apply a different map, a 
sequence of random phase shifts on the coin. We assume 

that at each step a random phase shift (9 a ) described 
by a unitary operator 



R(0M) = e 40 -" 5 Pa 



(3.1) 



a£D 



on the direction register is applied to a particle. This 
procedure, in effect, is equivalent to an application of 
another (random) coin on the whole direction register. 
Each random sequence E = {R(9^)}^L 1 generates a 
different quantum walk 

n 

U(9W) ■ ■ ■ U(8W) R{9^>)&j . (3.2) 

m— 1 

As we shall see later on, particular random walk proba- 
bility distributions associated with different sequences of 
random phase shifts do not differ significantly. 

What is significant on the QW-RPS is that it has mix- 
ing properties similar to the classical random walk (the 
linear growth of variance) and in particular for dimen- 
sion = 2, the mean probability distribution is exactly the 
same as for the classical random walk. Before we prove 
both statements, we derive the formula for the disper- 
sion of the probability distribution for QW-RPS, using a 
generalized Grover (i.e. permutation symmetric) coin. 

The generalized Grover coin is (cf. Eq. Q2.2|) l 

G r , t =r^P a + t J2 \a)(b\ = (r-t)l + t\D\P s . (3.3) 



a£D 



a.beD 

ajib 



where P a = \a)(a\. This operator is unitary iff the fol- 
lowing relations hold: 



M 2 + (\D\-i)\t\ 



1 



(\D\ -2)\t\ 2 + r*t + rt* = 



(3.4) 
(3.5) 



The QW-RPS may be thought of as a sequence of ran- 
dom operators U(ff) such that 



U r , t {0) 



E T a ® Pa 



aeD 



(1®C(0)) , 



and 



d(6) = 1® E el6a P* G r 



(3.6) 



(3.7) 



aeD 



Here 6 — (9%, . . . , 0\d\) is a sequence of independent ran- 
dom real variables. Hence we actually get a sequence 
of random operators {U(9^)}j which creates the QW- 
RPS. 



One step of a QW-RPS is given by: 

u r , t (dW) = s ■ cyW) = ( E T Q1 ®K)( ai 



aiED 



X [l ® E e W «\t^\D\\a)( S \ + (r- t)P a ) 

aeD 

= J2 T ^ ® e<> (VPIkM + (r - t)P ai ) 



ai£D 



(3.8) 



For the chain of n evolution operators of QW-RPS we 
obtain 

Lemma 3. Let {U r ,t(9^)} °°^ 1 be a sequence of random 
operators according to Eq. i^.b}) . Then 



ai ,...,a n £D 

® e 4( ^ )+ - +e - ) [(r - t)|oi) (a n \ + ty/\D\\ai) (s\ 

X H(oi, . . . , On)|oi) , 

where 

n-l 

S(oi, ....On) = Yl [ t+ ( r - t)8 aj ,a j+ i 
3=1 

Proof. By induction. 



(3.9) 



(3.10) 



□ 



It can be analyzed how a specific QW-RPS evolves 
given a specific initial state |^o) = |0) ® \s). 

Lemma 4. Let {U(9^)} C *L 1 be a sequence of random 
operators according to Eq. \3.b)) . Then 

\i>(Q,n)) = U{9^) ■ ■ - U{9^)\Q) ®\s) = 

= E |0 + e ai +--- + e a J«) 



ai,...,a n £D 



-lU) 



+•••+(?£ 



^3(oi, . . . , On) 



r - t 



\D\t 



Proof. From Lemma |3J 



\ai) . 
(3.11) 

□ 



The probability distribution of |^(0, x, n)) shall be de- 
rived by projecting it onto P x ® 1 and tracing over the 
coin Hilbert space. Hence by setting 



a = (oi, . . . , a n ) 
a' = (a'x, . . . , a' n ) 
9(a) = (9 ai + ■ ■ ■ + 9 a J 

we get 

P(9,x,n) = || (P x ® l)^(0,n)|| 2 



(3.12) 
(3.13) 
(3.14) 



r - t 



D\ 



\D\t 



E 



= i(0(a)-O(a')) 



3(a)S(a')*<aiK> , 



a,a'€i) n 
a=x=a' 



(3.15) 
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where 9(a) is the sum of the sequence of independent 
random variables for each a and S is the same as in Eq. 
(12.17(1 . The symbol a = x means that + S"=i e aj = x. 
It is clear that P(0,x, n) — P(S,x, n)* since we are 
summing over all tuples a, a'. 

Eq. 1(3.15(1 depends on 9, which is an event generating 
n\D\ random variables {9 aj : j = 1, . . . , n, a G a}. We 
can split Eq. 1(3.15(1 into two parts: 



P(6,x,n) = 



3(a)| 2 




+ 



(3.16) 

To obtain the mean probability distribution, we integrate 
over the random variable 9, assuming uniform distribu- 
tion for all phases of 9: 

i-2tt J\D\na 

(P(9,x,n)} e = / _ 7ii -P(e,x,n) 



o (2tt)I- d I«" 



r - 1 



E n 



-'a)| 2 



aeD n 
a=x=a' 



(3.17) 



The terms coming from ^2 a ^a' cancel out. Now the 



a x a 



mean probability distribution depends only on the term 
H(a). For d = 2 and Grover coin (i.e. r = rpr — 1, t = 

p), this term is the product of ±2, and |S(a)| 2 = 

2 -2(n-l)_ Xhen 



(P(e,x,n)) e = ^- E 1 



(3.18) 



; i x 



The sum over all paths in Eq. 1(3.18(1 is the same as the 
sum of all classical paths. The constant is the product of 
the probability to take any individual direction at each 
step. Hence the mean probability distribution of QW- 
RPS for dimension 2 is the same as the CRW in 2 dmen- 
sions. We easily show that the probability distribution 
resulting from Eq. (|3.18|) is normalized: 



E 



1 



<P(e,x,n)>e = -^ £ 1= I 

atD" 



(3.19) 



It is an interesting question whether the equivalence of 
QW-RPS with CRW in two dimensions is merely a co- 
incidence, or whether the same result applies for higher 
dimensions with a modified Grover coin. 

Now we arrive at the conclusion that the averaged 
probability distribution of QW-RPS with the generalized 
Grover coin G r ,t, is 



(P(9,x,n))e = 




)| 2 • (3-20) 



a£_D" 



This equation is one of the main results of our paper, as 
it shows that the probability distribution of a QW-RPS 
corresponds to a classical sum over paths. 

In order to study specific properties of the mean prob- 
ability distribution let us consider its dispersion that is 
defined as 



V(D,n)= ]T<P(e,x,n))e|x| : 



(3.21) 



The dispersion corresponding to the mean probability 
distribution is evaluated in Appendix A [see Eq. l(A15(l ] 
from which we can derive the following theorem: 

Theorem 5. The dispersion of QW-RPS with a gener- 
alized Grover coin with coefficients r,t for n > 2 is 



V(D,n) 



1 



l-|r| 2 + |£| 2 



(n-2) + (M 2 -|i| 2 ) 



(3.22) 



We want to find the coefficients r, t such that the above 
equation is identical to the probability distribution of a 
CRW. The sufficient condition is that all the terms in 5 
have the same absolute values. It is obvious that this is 
equivalent with the following requirement 

Theorem 6. The dispersion of QW-RPS with Grover 
coin G r t is identical to that of a CRW if and only if 

M = 1*1- 

Comparing this requirement with Eqs. 1(3. 4|l - ((3.5(1 it 
follows that the condition \r\ — \t\ is satisfied only in two 
cases: 



• dimension 
? + kir 



1, r 



73 c 



t 



V2 e ' 



a 







• dimension = 2, r = 



\ t ■ 



\e l! \ a- P = 7r+2/br 



The average probability distribution of QW-RPS with 
generalized Grover coefficients r, t is actually identical 
to the probability distribution of a CRW with memory 
(CRW-M). We define a CRW-M as the random walk of 
a particle on a <i-dimensional lattice, whose direction is 
changed at each step, depending on the direction from 
which it came. The CRW-M is given by the sequence 
{(x„, an)}™!, where x n 6 T, d is the position of the par- 
ticle, e an = x„ — x„_i is the unit vector in any of the 
2d directions, and ao is preset. One step of the CRW-M 
is given by the transformation (x„,a„) — > (x n+ i,a„ + i) 



such that x r . 



where 



Prob(a„+i|a„) = 



«n+l — «n 

otherwise 



(3.23) 



Beginning with xo = and ao in the uniform mixture 
of |D| = 2d directions, the probability Prob(x„ = y) is 
given by the sum of all paths from to y, each weighed 
by the product of terms |r| 2 , |t| 2 , depending on whether 
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the path continues in the same direction for two consec- 
utive steps. The sum of the amplitudes is the same as 

in Eq. (|3.20|l . weighted by the factor = ji, which 

corresponds to the mixture of different values of the ini- 
tial direction of the walker ao (see Eqs. (|3.4H . 113. 5(1 1. 

The memory effect in the QW-RPS with Grover coin 
is due to different absolute values of the coefficients r, t. 
Although Eq. (|3.2UI) shows that a QW-RPS using a sym- 
metric coin yields a probability distribution which corre- 
sponds to a CRW with memory (QRW-M) , it is also inter- 
esting to consider what QW-RPSs correspond to CRWs 
with no memory. Using the Fourier coin instead, and 
using random phase shifts will give the mean probability 
distribution equivalent to a CRW. 

The Fourier coin is defined by the operator of a d- 
dimensional Fourier transform 

^^EE^ww ( 3 - 24 ) 



j=0 k=0 



One step of d-dimcnsional QW-RPS with Fourier coin is 
defined by unitary operator 



U F (9) = 



LaeD 



(1<8F(0)), (3.25) 



with 



Obviously, 



a£D 



a\ . . ..,a n £D 
® (P ai F\D\ ■ ■ ■ Pa n F\D\) 



(3.26) 



(3.27) 



As before, we study the mean probability distribution 
induced by QW-RPS with Fourier coin and we conclude 
that the cross terms to vanish. It is straightforward to 
prove (cf. Lemma |2l : 



Lemma 7. 



P ai F w ...P an F w \s) 



°a n ,\D\ 



|D|(n-i)/a 



(3.28) 



x e 27ri(aia 2 H \-a„-ia n ) / \D\ i \ 



Now the probability distribution after n steps, with 
the initial state |0) <£> \s) and the sequence O of random 
phases is (cf. Eq. I3.I5|) 



Pp(6,x,n)= J] \Z F (*)f 



a.eD n 
a=x 

E 

a^a' 
a=x=a 



o t(fl(a)-0a')n 



5 F (a)H F (a')*(oi|a , 1 ) , 



(3.29) 



where 



"f(oi, ■ ■ ■ , ««) 



3 27ri(aia2H ha„-io„) /\D\ 



\D\i n - 1 )/ 2 

(3.30) 

By averaging over all sequences of random phases O (cf. 
Eq. (|3.f 7[) ) the second term in Eq. I|3.29|l vanishes and 
we get 



(F F (9,x,n))e = — ^ ^ S a n ,\D 



(3.31) 



Notice that despite the symmetric initial coin state |s), 
there is an asymmetry manifested in the term 5 a tj}\. 
This is due to the fact that after the initial coin toss, 
the coin register is in the state F\jj\\s) = \\D\). We can 
symmetrize the evolution by taking the initial coin state 
to be .fX-Js) = \\D\). Then we can compute that the 
probability distrubution has the same form as Eq. I|3.29|l . 
with Si? replaced by 



e 2Tri(a 1 a 2 A \-a n - 1 a n )i\ 



\ D \n/2 

yielding the mean probability distrubution 
(-PF, sy m(S, x, n))o = 1 • 



(3.32) 



(3.33) 



This is exactly the form for the probability distribution 
of a (memoryless) CRW, the sum being over all paths 
from to x, weighed by r^r. 



B. Numerical Results 

To complement our analytical results, we plot the dis- 
persion of the average probability distribution of the 
QW-RPS for the Grover coin: 



V Q (D,n) 



E 



x| 2 Tr(P x |^(e,n))(V-(e,n)|) 



(3.34) 

The results are shown in the following figures: in Fig. 1 
we plot the dispersion H3.34|l as a function of the number 
of steps for a two-dimensional (d = 2) QW, CRW and 
QW-RPS, respectively. 

Comparing the three corresponding lines we find that 
the dispersion in the QW grows quadratically with num- 
ber of steps 01 ( see solid/red line in Fig. 1). This is 
in a sharp contrast to a classical random walk for which 
the dispersion is a linear function of the number of steps 
(see dotted/blue line in Fig. 1). The lines for QW-RPS 
and CRW overlap. The dispersion of either grows lin- 
early with the number of steps n (see dashed/blue line in 
Fig. 1). 

In Fig. 2 we plot the dispersion l|3.21|l as a function 
of number of steps for a three-dimensional (d = 3) QW, 
CRW and QW-RPS, respectively. 
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FIG. 1: (Color online) Dispersions of the probability distribu- 
tions corresponding to the n steps of quantum walk (solid/red 
line), the (memoryless) classical random walk with equal 
probabilities of step in any direction (dotted/blue line), and 
the quantum walk with random phase shifts (dashed/green 
line) for a two-dimensional system. The initial state of the 
quantum system system is described by a vector |?/>o} = 
|0) ® \s) and we assume the Grover coin Cg- The quanti- 
ties for QW-RPS were obtained by generating 50 evolutions 
of QW-RPS with respective dispersions of probability distri- 
butions, and by averaging over them. 



FIG. 2: (Color online) Dispersions of probability distributions 
corresponding to the n steps of quantum walk (solid / red line) , 
the (memoryless) classical random walk with equal probabili- 
ties of step in any direction (dotted/blue line), and the quan- 
tum walk with random phase shifts (dashed/green line) for a 
three-dimensional system. The initial state of the quantum 
system system is described by a vector \if)o) = |0) ® |s) and 
we assume the Grover coin Cg- The quantities for QW-RPS 
were obtained by generating 50 evolutions of QW-RPS with 
respective dispersions of probability distributions, and by av- 
eraging over them. 



As in the two-dimensional case, the dispersion of the 
probability distribution of the QW grows quadratically. 
The dispersion of the classical random walk is a linear 
function of number of steps and it does not depend on 
the dimension of the random walk. Interestingly enough, 
the dispersion of the quantum walk with random phase 
shifts is again a linear function, but unlike in the two- 
dimensional case, for d > 3 the linear growth of the dis- 
persion is faster than in the classical case. 

The same conclusions can be derived from our simu- 
lations of quantum walks in four-dimensional space (see 
Fig. 3). 

In Fig. 4 we plot dispersion of probability distributions 
for quantum walks with random phase shifts as a number 
of steps for various dimensions d — 2,3, 4. We generated 
50 evolutions of QW-RPS for each dimension, and av- 
eraged over the respective dispersion generated by each 
evolution. 

We can conclude, that as the dimension increases, the 
linear growth of the dispersion also increases. 

We have shown that the introduction of random phase 
shifts causes the transition of a QW to a (quasi-)classical 
random walk. In our previous discussion we have consid- 
ered random phases to be uniformly distributed in the in- 
terval [— 7r, 7t]. Here we briefly analyze a situation when 
phases of random shifts are distributed according to a 
normal distribution that is peaked around the phase zero 
and with the dispersion a. When the dispersion is zero, 
(i.e. a — 0) we recover the QW (see Fig. 5), while for 



large a, we obtain uniform distribution on the interval 
[— 7T, 7t] and the CRW is obtained. The results are shown 
in Fig. 5. This analysis clearly shows the quantum-to- 
classical transition for quantum walks which is generated 
by random phase shifts. As the phase shifts become more 
random, the walk becomes more classical. 

IV. CONCLUSION 

We have shown that by shifting the amplitudes of the 
coin register in a quantum walk by random phases, we 
can obtain the classical behavior of the quantum walk. 
For a Grover coin, the mean probability distribution of 
such a walk is equivalent to the CRW with memory and 
a biased coin; for the Fourier coin, the mean probabil- 
ity distribution is equivalent to the memoryless CRW 
with an unbiased coin (given an unsymmetric initial coin 
state). 

The results underlying Fig. 5 also show how the transi- 
tion from QW-RPS to CRW occurs when we increase the 
dispersion of the normal distribution of random phases 
(for the Grover coin). Our results are in a way comple- 
mentary to a standard quantization procedure in physics. 
Specifically, classical dynamics of physical systems can be 
canonically quantized, so it is clear what is the quantum 
version of a classical process. On the other hand, quan- 
tum walk is not obtained by a canonical quantization 
procedure from a classical random walk. It is simply de- 
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FIG. 3: (Color online) Dispersions of probability distributions 
corresponding to n steps of quantum walk (solid/red line), the 
(memoryless) classical random walk with equal probability of 
step in any direction (dotted/blue line), and the quantum 
walk with random phase shifts (T>o(D,n): dashed/green line) 
for a four-dimensional system. The initial state of the quan- 
tum system is described by a vector \ipo) = |0) ® |s). We 
assume the Grover coin Cg -The quantities for QW-RPS were 
obtained by generating 50 evolutions of QW-RPS with respec- 
tive dispersions of probability distributions, and by averaging 
over them. 



FIG. 4: (Color online) Dispersion of QW-RPS processes for 
different dimensions for n steps. We see that these dispersions 
T>o(D,n) are linear functions with gradients that depend on 
the dimensionality of the system under consideration. Only 
the case d — 2 coincides with the classical random walk. The 
quantities for QW-RPS were obtained by generating 50 evo- 
lutions of QW-RPS with respective dispersions of probability 
distributions, and by averaging over them. 



fined by a set of instructions that govern the evolution of 
the quantum walk. Therefore it is of importance to know 
what is the underlying classical process. This underly- 
ing process can be reconstructed either by measuring the 
coin at each step (cf. ;6j), or when the quantum walk 
is subject to random phase shifts that totally suppress 
quantum interference between different evolution paths. 
As a result of the suppression of the quantum interfer- 
ence, the classical random walk that corresponds to the 
underlying quantum walk emerges. 
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FIG. 5: (Color online) Dispersion T>o(D,n) of n steps of QW- 
RPS in dimension 2 with random phases normally distributed 
around with respective standard deviation. We see a "con- 
tinuous" transition between the QW and CRW as function of 
the standard deviation of the random phase distribution. 



APPENDIX A: DISPERSION OF QW-RPS WITH 
GROVER COIN 



The dispersion reads 



Starting with Eq. l|3.20f> we can evaluate the dispersion 
V{D,n) of QW-RPS with generalized Grover coin G r .t- 



V(D,n)=Kj2 |e ai +--- + e a J 2 |S(a)| 2 (Al) 



aED" 



9 



with K = 



>\D\ 



\D\t 



Turning Eq. (|A1(1 into the 



recursive relation we obtain 
V{D,n)=K Y, {l e -! 



+ 2(e 0l +--- + e 0n )-e a „ +1 +l} 
x |S(oi, • ■ • ,a„)| 2 |S(a„,a„+i)| 2 . 



(A2) 



We find that E a „ +1 eD a„+i)| 2 = V\ 2 + (\D\ - 

l)\t\ 2 = 1 for all a n € D. Hence the first and the 



last terms in the braces contribute to Eq. (|A2f) with 
D(D,n) + 1. The middle term has the form 

2K Y |S(a)| 2 (e 0l +---+e a „)- ^ |S(a„, a„+i)| 2 e a „ +1 

a.ED™ a n + 1 eD 

(A3) 

In the sum over a n +i in Eq. (|A3|) . we can keep just the 



terms a„+i such that e Qll+1 is parallel with e ai +• ■ - + e aii . 
The remaining a n +i-S cancel out, since for each such a n +i 
there is a^ +1 such that e Qra+1 + e' a — 0. Hence the 
second term of Eq. I|A3[1 can be rewritten as 



2K{\r\ 2 -\t\ 2 ) ]T |H(a)| 2 (e 0l +--- + e J-e a „ 



(A4) 



2K(\r\ 2 - \t\ 2 )R n . 



The expression for V{D,n) reads as 

T>(D, n + l)= V(D, n) + l + 2K{\r\ 2 - \t\ 2 )R n (A5) 
with 

|2 

(A6) 



V(D,2)=K J2 \^{a l ,a 2 )\ 2 \a l +a 2 \ 

= X{4| J D|H 2 + 2(|^| 2 -|^|-l)|t| 2 } 



The expression R n can be rewritten into recursive equa- 
tion: 



R n +i = (\i 



with the initial condition 



\t\ 2 )Rr, 



1 

K 



R2= Y l"( a i ; a 2)| 2 (ea 1 +e Q2 ) -e Q2 



(A7) 



(A8) 



Eqs. I|A7(1 and (|A8|) can be solved to obtain 

(|r| 2 - \t\ 2 ) n - 2 {KR 2 \r\ 2 - KR 2 \t\ 2 - KR 2 + 1) - 1 



Rn — 



K{\r\*-\t\*-l) 



(A9) 



Solving Eq. l|A5p and collecting the terms from Eqs. 
(|A"6)l . jSHt , and fXojt we obtain 

V(D,n) 



where 



(|rp-|tp)(l-|r| 2 + |t| 2 ) 2 
{(n - 2)$ - 2(|r| 4 + |<| 4 ) + 4|r| 2 |i| 2 + (A10) 

(2 + 7 ? )(|r| 2 -|t| 2 )"+r ? [|t| 2 -|r| 2 ]}, 



C=M 2 -|r| 6 



and 



rj = 2 



i| 2 + 3|r| 4 |t| a -3|r| 2 |t| 4 + |t| 6/ (All) 
r + (|D|-l)f 



\D\ 



(M - |*| - 1) 



{2|i%| 2 + (|D| a -p|-l)|i| 9 } 



(A12) 



We may assume that r is real and t = |i|e M . Solving 
Eqs. I|3.4[) and l|3.5(l we obtain 



t = 



1-r 2 



1/2 



a = ± arccos 



2r 



(2 -pi 



1-r 2 



2 \ 1/2' 



(A13) 
(A14) 



where |D| > 4, < r < 1. Obviously, < |r|, |i| < 1 



and Irl 2 - lil 



PI 



rji^j-. Eq. I|A10|) contains only two 

terms dependent on n: (n — 2)£ and (|r| 2 — |t| 2 )". The 
latter goes to as n — > oo, hence we get (n > 2): 

1 + Irl 2 -It' 2 



V(D,n) 



1-H : 



^(n-2) + 0((|r| 2 -|t| 2 r) . 

(A15) 



= 2IDIIH- 



|D| 2 -|£|-l)k| 2 
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